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I. 



INTRODUCTION 



One important design consideration of any naval vessel is 
the habitability of the craft for all possible environmental 
conditions encountered in its operation. For the CAB type surface 
effect ship, it is important that the vibrations induced by the 
high speed operations over the full spectrum of the sea surface 
wave profile do not produce a great discomfort factor for the 
crew. In particular, the possibility of large vertical accelerations 
for prolonged periods of operation should be avoided. 

The primary purpose of this study is to obtain data relating 
the effects of plenum air flow leakage, and air influx rates on 
the vertical motions of The Captured Air Bubble craft. To do 
this several models and several kinds of "simulation tests are 
used. The models used are a Six Degree of Freedom computer 
model (originally written by Oceanics, Inc., and modified for the 
XR-3 craft at USNPGS) , and the simplified linear model for heave 
motions only (written at USNPGS, see Report NPS62-77-002 , Ref. 2). 

The simulation tests included sinusoidal sea frequency 
response tests, complex sea frequency response tests, and step 
weight removal transient tests. 
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II. SIJyLPLIFIED XR-3 MODEL 



In order to investigate the characteristics of the vertical 
motion of the CAB craft, it is necessary to obtain a mathematical 
model that can readily be used for analytical purposes. This was 
the principal objective in Reference 1 in which a simplified 
heave-only model of the XR-3 was developed. For the convenience 
of ready reference, the development of this model has been 
repeated in Appendix A. 

This simplified XR-3 model has proven to be quite useful in 
developing a scaling method for the heave motion equations as 
reported in Reference 2. The linear model used was described 
by the state variable vector 

X = (62 , hz , < 5 V^) 

and the output variable , 

y = <sPb 

as given in Appendix A where 

xi = 6 z, incremental heave (draft) 
xz = 6 z, incremental heave rate 

X 3 = <SV^, incremental air mass volume with 6 V^ = 

y = <SPj^, incremental plenum pressure 

The step weight transient used both in Reference 1 and Reference 2 
had the state equation 
0 

X = A X 
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with initial conditions 



x(0) = [6z(0) 0 6V^(0)]’^ 

which was brought about by the sudden removal of 10% of the craft 
weight. The characteristic equation for the system 

det(sl - A) = 0 

is shown in Appendix A to be of third order. In this report, the 
sinusoidal sea input transfer function for this model is developed 
in the following section. The step weight transient is investigated 
in Section IV to correlate the results of the sinusoidal sea input 
investigation . 
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III. SINUSOIDAL SEA INPUT 



The characteristics of the fully developed sea surface can be 
approximated by a long crested wave obtained by the addition of a 
finite number of sinusoidal waves of various amplitudes as described 
in References 3 and 4 and siimmarized below. 

For the CAB-SES frequency response analysis , it is necessary to 
develop a transfer function for the sinusoidal sea input. What 
follows is a development of the transfer function for the heave motion 
equation of the simplified XR-3 craft. 

A Development of the Simplified XR-3 Transfer Functions 
1. Surface Wave -Profile 

The surface wave profile for both deep and shallow waves can be 
represented by a regular swell or wave: 

WAV(x,t) = A sin[k (x + V t) ] 

n w 

where: A = surface wave amplitude (half-height 

from crest to trough) (ft) 

= the "wave number", 2 it/L^ (ft~M 

L = wave length (ft) 

w 

= wave velocity or "celerity" (ft/sec) 

^or deep water, = 

A more convenient form can be obtained by using circular (inci- 
dent) frequency, = 2 it/T^ where T^ is the period, that is, the 
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time required for the wave to travel one wave length. From the 

/2-rrL 

equation for ^ 

then 

w \ w / 

thus , 



WAV(x,t) = A sin(k x + (O.t) 

n 1 

When observed at a fixed point in space with x = 0 



WAV(0,t) = A sin(w^t) 

Alternately, if the wave profile is observed at t = 0 , 



WAV(x,0) = A sin(k x) 

n 



The slope of the wave surface is obtained by differentiation. 



|j[wAV(x,0)] = k„. 



A cos (k x) 
n 



and the maximum slope occurs when k^x 



= 0 , IT , etc. , thus 



- , TTh 

max slope = k A = :r^ = :p— ^ 
n L ^ 



w 



w 



where = length of the wave 



h„ = 2A = wave height from hollow to crest, 
w 

2. Wave Input for the Heave Equation 

In the case of the CAB-SES, the wave introduces (1) a change in 
the plenum volume and thus a change of plenum pressure and (2) a 
change in sidewall immersed volume and thus a change in the forces 
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of buoyancy. The additional forces due to the seals will be neglectec 
in this study since they largely affect pitch motion which is not 
considered in this model. Since the objective is to obtain transfer 
function of the system to be analyzed for the condition of a forcing 
function that is sinusoidal in time, it is necessary to recast the 
wave profile equation from a function of displacement and time to a 
function of time only. This is readily possible for the heave-only 
model since all the forces along the hull that result due to the waves 
can be lumped together into a single force acting at the center of 
gravity (CG) which is the same point as the center of pressure (C.P.) 
for this simplified XR-3 model. 

Given the wave profile 

WAV(x,t) = A sin(k x + co.t) 

n 1 

and the assTomption that the craft is heading into the wave, the wave 
height, WAV(x,t), as seen by the ship hull is at the apparent or en- 
counter frequency, + ^—) where V is the ship velocity and 

w 

V is the velocity of the wave, thus 
w 

WAV(x,t) = A sin(k x + <o t) 

n e 

= A sin / 1 + CO t\ 

® ; 

w ' 

Use of the identity for the sine of the-sum-of-two-angles yields 

[ CO . X CO . xT 

sin(co t)cos + cos ( CO t) sin 

e V e V I 

w w J 

The differential volime of the plenum due to the wave input is 
d[5Vp(x,t)] = (Wp) [WAV (x, t) ] dx 
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where 



5V is the incremental plenum volume due to wave. 
P 

W is the width of the plenum chamber. 

P 

dx is the differential length of the plenum. 



The total incremental volume of the plenum due to the wave 
can be established at any time t by integration over the plenum 




It can be observed that the first term on the right hand side 
of the equation when the limits are applied will yield a double- 
amplitude sinusoid (since it is an odd function) and the second 
term will yield zero (since it is an even function) . Thus, 



%(t) = 



V ' 
^ p w 



sin 




w /J 



A sin(o) t) 
e 



9 



Now, for each regular wave input with a wave length of and 



2irV 



an incident (circular) frequency of — — where = 



w 




w 



the 



coefficient in the brackets above is constant and thus the incre- 
mental variation of plenum volume is the result of the sinusoidal 
input, 6W = A sin(oj^t) . 

The transfer function of the incremental plenum volume, 5V 



to wave height, 6W becomes. 



5V 



«w ' = 




where the coefficient ki (co^) is a function of the incident frequency. 

Following the same path of development, the transfer function of 

each sidewall volume, 5V , to wave height, 6W is, 

s 




ko (Wj, ) = 




where W = width of each sidewall . 
s 

The transfer function of the total incremental sidewall buoy- 
ancy force, 6F to wave height, 6W is 
s 



5F 

c 

W 



= ko ) (2p^g) 



Then 



5F 



sm 



ko (oj^) 



2p g 
M 



where 6F = 6F /M 
sm s^ 
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3. System Block Diagram and Transfer Function 

The linear heave motion equations were developed in Ref. 1 and 
are repeated here for convenience in Appendix A. The block diagram 
for this linear system is shown in Figure 1 where the output is the 
incremental variation in heave (or draft) 5z, and the input is the 
sinusoidal variation in wave height about the calm water level of 
draft. Note that a positive value of wave height reduces the effec- 
tive pleniam volume, , and therefore requires a negative sign in 

the signal transmission as shown in Figure 1. 

Also it should be noted that a positive value of wave height 
increases the wetted sidewall volume, , which cause a larger 

upward force and therefore requires a negative sign in the signal 

transmission. The effect of added mass, 5 f , was also included in 

am 

this model. However, in the development of the transfer function of 
C.G. acceleration to wave height as well as plenum pressure to wave 
height, it was shown that the effect of added mass can be neglected, 
as a reasonable first approximation and therefore added mass effects 
were not included in the analysis that follows . 

The algebraic computation of the wave input transfer functions 
can be performed either by block diagram manipulation or by applica- 
tion of the signal flow graph method. The latter approach was used 
with the details contained in Appendix B. The resulting transfer 
functions are of the following forms 



Sz 

(SW 



- (a s + b) 

s^+cs^+ds+e 



(Draft/Wave height) (1) 
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FIGURE 1. BLOCK DIAGRAM OF THE LINEAR HEAVE EQUATION 

FOR SINUSOIDAL WAVE INPUT 




« • 

6z 

6W 



-s^ (a s + b) 

+ cs^ + ds + e 



'C.G. acceleration \ , - > 

Wave Height ) 




(a ' s + b ' ) 

+ cs^ + ds + e 



/ Plenum Pressure 
\Wave Height 



The coefficients a,b,a' and b' are functions of the incident 
wave frequency since they contain the ko and ki terms of the plemam 
volume-wave height transfer function and sidewall volume-wave height 
transfer function developed in secion III-A-2. It is possible to 
simplify this functional relationship in ko and ki by a convenient 
approximation which is carried out below. 

From Appendix B, 



a = ko + kik2ksk4 



where ko = 



ki = 




('2W V \ /o) .L \ 2p g 

s w\ . / 1 p \ 0 )’ 

‘ 2 ^) -JT 

''oj . L 



sin 






2V 



w 



k2 = 1/Vj^{0) 
k3 = YPj^lO) 
k4 = Aj^/M 



For the XR-3 craft operating at 
quency, , can be approximated, as 
frequency range of interest as, = 
g/oi^ , it is possible to apply these 
to yield 



30 knots , the encounter fre- 
shown in Appendix C, over the 
2(0^^. Recognizing that = 
substitutions into ko and ki 
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ko = 



2A p g \ 

s \ sin ao3 



M 



aco 



and 



ki = 






sin aoi 
ao3 



where 



A, = W L 
b p p 



a = L /4g 
P ^ 



p s p 



CO = 0), 



Observe now that ko and ki are now composed of constant terms 
multiplied by the sine function A plot of this function 

over the encounter frequency range of interest is given in Figure 
2. It should be noted that the magnitude of this function is very 
nearly unity up to 2 rps then attenuates in a manner that approxi- 
mates a single pole at 8.4 radians/sec up to 15 rps after which it 
drops at a faster rate and then a lobe of increased magnitude 



occurs with a peak at approximately 30 rps and repeated attenuated 
lobes at multiples of the plenum length as reported in reference 7. 



The coefficient 



a = 



2AjP„g ^ YPb<“>'^b 



M 



V, ( 0 ) M 
b 



and in like manner 
b = 



2AsPwg YPj3(0)ks 



M 



a' = b' = 






Sin aco 
ao3 



sin aco 
ao3 



. ''b<“' 



sin ao) 
aco 
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RANGE OF u VALUES FOR WHICH THE HEAVE GAIN APPROXIMATION ISVAUO 
MAXIMUM RANGE *4 



a 




where the added mass coefficient, kg , that originally appeared in 
coefficient b has been set equal to zero. 

The transfer functions now become 



where 



and 



6z _ -a(s + b/a) 
6W 3 ^ 2 



(lb) 



s + cs + ds + e 



6 z _ -as (s + b/a) 

s^ + cs^ + ds + e 



(2b) 



'^^b _ a* s (s + 1) 

<SW 3 2 

s +cs^ +ds + e 



(3b) 



a = 



a = 



b = 



M 



yPv. (0 ) A, V 
‘ b b \ sin aoj 



M / aco 

M V„(0) 
m 



sin goj 
ao3 



(since >> 2A p g) 

' b b 



sin g(^ 
ao3 



a * = 



V^b"” 

v^,(0) 



sin ao) 

ao3 



ke = I n k + 
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V^(0) 
m 

“ M V, (0)M 

yP ( 0 ) A ^ 

■ v^roTi T 2*ppw9> 

Vb(0) 

YP^(0)ks f2A^P„g\ 

« - V^(0) \ M ) 



The value of b/a is dependent upon the magnitude of the air 
flow rate term, ks , thus 

. 2A p gyP, (0)k „ 

b _ s^w^ ' b 6 M 

a MVj^(O) * (2A^p^g + yPj^(O)Aj^) 



with added-mass coefficient 
kg = zero 



which will reduce by the above approximation to yield, 



b 

a 






kj 



Using the parameter values given in Section IV, the transfer 
functions are 



6z 

5W 



■-3200. (s + .045 ks) 
s^ + 11.3 ksS^ + 2441.8s + 145.3 



sin (.1550)) 
(. 1550)) 



5z _ 2 <Sz 

6W ® 6W 




2284 s^ (s + 1) 

+ 11.3 ks + 2441.8s + 145.3 



sin (. 15503) 
(1.55o)) 



17 



where 



a = 



b = 



a' = 



c = 



d = 



2 (^) 2 (32.17) + 1.4 (2140.8) 200 



188.06 



'"2(j^) 2(32.17) 1.4(2140.8) kg' 



sin ( 



20 



4 ■ 22 . 11 ‘ 



0 ) 



(5/32.17) 03 



188.06 (264.63) 
200(1.4) (2140.8) I sin aco 



sin ao) 
aco 



262.44 



.8) j 



aco 



1.4(2140.8) 1 
264.63 J 



2 (|=) 2(32.17) ^ 



188.06 



188.06 



e = (11.326) (12.83) kg 



Using the value of kg = 6.28, the transfer f;mctions become 



sin ( . 155co) 

6z _ -3200 (s + .28) ( .155co) 

6W (s + .38) (s + 35.1 + j 34.4) (s + 35.1 - j 34.4) 



6 z _ 2 6z 

6W ® ^ 



rp sin(.155 ) 

b ^ 2284 s^ (s 4- 1) (.155 ) 

6W (s + .38) (s + 35.1 + j 34.4) (s + 35.1 - j 34.4) 



where the factored form of the denominator was obtained from Section 
IV A. 

Figure 3 shows the magnitude spectrum for the three transfer 
functions of the simplified XR-3 model for the value of kg = 6.28. 
The curve for acceleration was plotted in g units and the pressure 
curve in incremental (gage) pressure values. 
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FIG 3. FREQUENCY RESPONSE OF THE LINEAR - 

SIMPLIFIED MODEL LOW LEAKAGE, K|=6.28 
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B. Calculation of the 6 DOF Model Frequency Response. 

Using the 6 DOF equations for the XR-3 as given in Ref. 5, 
the output response to a small amplitude (10% of the average draft) 
sinusoidal head-on sea was obtained on the IBM 360/67 computer. 

Two types of wave inputs were used; a single frequency input 
(regular sea) and a complex frequency input (complex sea) . The 
band of encounter frequencies was extended beyond that used in 
Ref. 6 to include the range 

1 < CO < 100 rps 
1. Regular Sea 

The procedure used in the computation of output response to 
a regular sea input was to compute 35 seconds with the adjustable 
thrust option to arrive at a steady state sinusoidal response. 

The average value of thrust, draft, pitch and pressure were then 
used as input for the operating condition of the final run. The 
steady-state final run was obtained with constant thrust at a 
steady-state speed of 30 knots. Data was taken only after at 
least 30 seconds of compute time to insure true steady-state 
conditions. The integration method used was a f ixed-s tepsize 
Runge Kutta which yielded about 50% saving in compute time over 
the adjustable R-K method with no observable loss in accuracy. 

The ratio of draft to wave height amplitude as well as 

C. G. acceleration and plemam pressure to wave height were computed 

for the 6 DOF equations with the results plotted in Figure 4 . All 
3 frequency response curves clearly show the function 



20 




-h H 

O 3AVM O 

S 3dflSS3ad 



-f- 

S 



o 



H 1 1 — 

o 1H9QH 3AVM o lO 1H9I3H 3AVM 

NOIlVd313XV9 0 - O liVdQ 



2 ./ 



lOQ 



lobes occuring at approximately 28 rps and above. A comparison 
of Figure 4 and Figure 3 clearly shows that the simplified XR-3 
is in good agreement with the 6 DOF results . For frequencies 
below 28 rps, the transfer functions can be approximated from 
the data of Figure 4 using the straight line approximation methods 
of Bode as follows : 

For draft, the frequency response follows very closely the 
curve of Figure 2 for the sine function thus 



(0 = oje 

For pressure, the response is approximated by the transfer 
function 



For C.G. acceleration, the response is approximated by the 
transfer function (in g units) as. 



These above approximations neglected the small amplitude 
variations that occur between 3 and 6 radians per second that are 
most likely introduced by the pitch mode coupling which is not 
accounted for in the simplified heave-onlv model of this study. 

2 . Complex Sea 

The complex or irregular sea conditions were approximated by 
the summation of 10 regular sea components as provided for in the 



5 z ^ sin OLOi 
5W ~ ao) 



where a = Lp/4g 
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6 DOF wave subroutine but with sufficient amplitude in each com- 
ponent to introduce an output response without exceeding the total 
height limitations for maintaining somewhat linear operation. The 
transfer functions amplitude ratios were computed at these dis- 
crete frequencies using the Fast Fourier Transform of Reference 4 
and the values plotted in Figure 4 as square points. It can be 
seen that each of the discrete frequencies produces values that 
are in reasonable agreement with the regular sea input results. This 
also has been confirmed by Booth in Reference 4. It should be 
noted that the complex sea input frequency response computation is 
a single run technique but limited by the following 

1) Discrete frequencies must by necessity be limited 
in number and so chosen that they are harmonics of 
the fundamental frequency. 

2) Each frequency must introduce a reasonable ampli- 
tude to produce a response which in total must be 
limited to maxim\am composite wave height. 

3) Steady-state speed must exist since encounter fre- 
quency must be a multiple of the fundamental fre- 
quency to reduce the effect of leakage (a problem 
with the Discrete Fourier Transform method) . 
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IV. PLENUM AIR FLOW RATE PARAMETER STUDY 



During the study conducted on scaling the heave motion 
equation as reported in Reference 2 , it was observed that 
the air flow influx as well as the leakage rate determined 
the characteristic response of heave motion for a step weight 
disturbance. It is of interest, therefore, to study the 
sensitivity of the craft dynamics to these parameters (air flow 
influx and leakage rate) . The characteristic equation of 
our simplified linear model is a function of the parameters, so 
a root locus study permits analysis of their effect on root 
locations. Alternately the transfer functions may be recomputed 
and the frequency responses plotted and analyzed. Finally the 
six degree of freedom simulation program may be exercised to 
provide both frequency response data and step weight removal 
transients . 

A. Characteristic Equation Root Locus 

The linearized equations (see Appendix A) yield the following 
characteristic equation 

s^ +cs^ +ds+e=0 

which in terms of the system parameters becomes. 
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represent the adjustable 



Now let = 

parameter for the characteristic equation which can be re-organized 
into the root locus form for a 7r-type locus , 





K 



(1.4) (2140.8) 
I 264.63 



2 ^ 2(1.99) (32.17) (75) (32.17) 



s^ + 



(4) (6050) 



-1 = 



1 



s<s^ + 



2 (75) (32.17) (1.99) (32.17) (1.4 ) (2 0 0) ^ (2140.8) (32.17) 



4(6050) 



(262.44) 6050 



Then the equation for the root locus becomes 



(11.33) (s^ + 12.83) 
s(s^ + 2441.8) 



1 



To carry out the root locus study numerical values are 
required for the parameter. We obtained the values used in 
the following study from the six degree-of-freedom program. 
These values were also used in the simulation studies for 
the 10% weight removed transient; thus the root locations 
can be correlated with the transient response. 
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An IBM 360/67 subroutine called ROOTLO was used to compute 

and plot the locus of roots with the results shown in Fig. 5 and 

in Table I. The values of K used in this study covered the 

x» 

range 



6 < < 22 

with the root locus of Fig. 5 covering the range 

5 < < 50 

From Fig. 5, and Table I, it can be seen that the roots of 
the characteristic equation for small values of are one real 

pole near the origin and a complex pair far removed from the origin. 
For larger values of the complex pair become real roots with 

one root approaching the root near the origin and the other moving 
towards negative infinity. Note in particular that the real root 
near the origin moves a relatively short distance during the one 
decade change of the parameter K^. 

In summary, it is noted that the real root near the origin 
is not very sensitive to changes in this range of values, 
whereas the two other roots are quite sensitive to changes in K^. 

Since K. is a function of air flow rate influx (nk_) and air flow 

Fi. 

rate leakage -/■= 7 -^ , it is clear that the oscillatory 

Pa V2Pb(0) 

modes in the ship dynamics depend on these flow rates. 
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Table I 



- Root Location for Various Values for the Characteristic 
Equation, (s+a) (s+S) (s+y)=0 







a 






8 








Y 


5.00 


0 




3 


49.4 


0 




3 


49.4 


0.0 


6,28 


- 35.1 


+ 


j 


34.4 


- 35.1 


_ 


j 


34.4 


- 0.38 


7.88 


- 44.0 


+ 


3 


21.5 


- 44.0 


__ 


3 


21.5 


- 0.47 


9.88 


- 81.3 




3 


0.0 


- 29.2 


__ 


j 


0.0 


- 0.60 


12.4 


- 119 . 


+ 


3 


0.0 


- 19.6 


_ 


3 


0.0 


- 0.77 


15.6 


-160 . 




3 


0.0 


- 14.2 


__ 


3 


0.0 


- 0.99 


19.5 


- 208 . 


+ 


3 


0.0 


- 10.4 




j 


0.0 


- 1.31 


24.5 


-267 


+ 


3 


0.0 


- 7.3 




3 


0.0 


- 1.82 


28.5 


-313 


+ 


3 


0.0 


- 5.3 


__ 


3 


0.0 


- 2.50 


33.2 


-366 




3 


0.0 


- 3.31 




3 


1.44 


- 3.31 


- i 1.44 


38.6 


-428 


+ 


j 


0.0 


- 2.84 




3 


2.23 


- 2.84 


“ i 2.23 


44.9 


-500 


+ 




0.0 


- 2.43 




3 


2.66 


- 2.43 


- j 2.66 
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B. Sinusoidal Sea Input Frequency Response 

Both the simplified linear model and the six degree of freedom 
program were used in this part of the study. Figure 6 shows the 
frequency response of the linear model for high leakage conditions. 
Comparison of Figure 6 with Figure 3 shows one effect of the 
leakage conditions: with high leakage the higher frequency 

motions are substantially attenuated. Note that for encounter 
frequencies less than 8 rad/sec (w^ < 8) the craft's vertical 
motion is not at all sensitive to air flow and leakage rates , 
but for w^ > 10 the CG acceleration, pressure and heave 
displacement are all substantially reduced for the case of increased 
leakage . 

Figxire 7 shows the frequency response of the six degree of 

freedom model for high leakage . Comparison of Figure 7 with Figure 4 

% 

shows the same general characteristics observed with the linear 
model: the response to frequencies above 10 radians per second 

are significantly attenuated when the air flow is increased. 
Frequencies below 10 radians per second produce somewhat larger 
responses when the airflow is increased. We believe this is a 
result of pitch motions, as discussed in the following paragraphs. 

C. Comparison of the Linear Model with the 6DOF Model Using the 
Frequency Response Plots 

One purpose of this study is to assess the value of the linear 
model as a tool for studying the CAB and predicting the salient 
features in its behavior. Alternately, comparison of results 
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between the linear model and the 6 DOF model (or experimental data 
when available) should lead to an understanding of phenomena 
which are peculiar to the CAB. On Figure 8 the c.g. acceleration 
curve of the simplified model and that of the six degree of 
freedom model are plotted, using the data for = 19.99. 

At low frequencies the six degree of freedom curve indicates 
much larger accelerations than the linear model. At high 
frequencies both models are in good agreement (the null frequencies 
are not exactly the same, indicating that the effective plenum 
length is not precisely the geometric length) . In the low 
frequency range we note that the six degree of freedom frequency 
response has a slight "dip" or anti-resonance for 4 < w < 6, 
while the linear model does not. This frequency range includes the 
natural frequency of the pitch mode, which is included in the six 
degree of freedom model, but not in the linear model. We know 
that the sinusoidal sea will excite the pitch mode, so we are 
led to a plausible explanation for the discrepancy between the 
results at low frequencies. It is simply this - with sinusoidal 
excitation in the range of frequencies about pitch resonance, 
the pitching motion will have large enough amplitude to increase 
the air leakage. With the bubble volume thus reduced the 
incompressible wave will cause higher air pressures and higher 
heave acceleration. At high frequencies, above the natural 
frequency in pitch, leakage is not increased and the linear model 
is a reasonably good one. 
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FIGURE 8 COMPARISON OF FREQUENCY RESPONSES IN ACCELERATION 
Ki =19.99 
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We also note that at low frequencies the slopes of the 
curves on Figure 8 are not the same. This can be attributed 
to variation in the leakage as a function of frequency. 

D. Simulation of Transient Response due to Step Removal of Weight 
The removal of weight by a step change is equivalent to a 
negative step change in acceleration. The transfer function for 
draft pressure and C.G. acceleration for the step weight removal 
can readily be found by Mason's Gain Rule using the signal flow 
graph given in Figure B-1 of Appendix B. The transfer function 
for draft becomes 




5a D(s) 



where 6z = incremental draft 

5a = incremental weight (acceleration) 



D(s) = s^+ 






K, s+ 



M 



YPb«»A^’ 

MV|^(0) 



S + 



^YPi3(0)' 



M 



For pressure. 



!!b .\V£LJ! 

5a D(s) 



And for C.G. acceleration where 5z = incremental acceleration: 



5a 



s^ (s 



YPb(O) 

D(S) 
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These transfer functions are for the linear simplified 
model, for which we have the root locations (factors of D(s)) in 
Table I with corresponding root locus on Figure 5. The transfer 
functions also indicate the presence of zeros , some of which 
are functions of the leakage, K^. 

In studying the weight removal behavior, simulation was done 
with the six degree of freedom model, and the results are 
conveniently analyzed if we consider Table I and Figure 5, 
remembering also that the simulation model contains additional 
pitch mode effects. Figure 9 is a set of 16 curves covering the 
range of air flow conditions from 6.1 < <21.99, but recorded 

for less than a 1.0 second duration after removal of the weight. 

This expands the scale of the initial transient, and permits 
correlation with Table I and Figure 5. Figure 10 is a set of 
12 curves covering the same air flow conditions, but recorded for 
a 5.0 second duration, which permits us to recognize the effect 
of pitch coupling. 

Consider first the c.g. acceleration curves of Figure 9 
(Al, A2, A3, A4) , noting from Table I that for each of these 
figures the value results in a small negative real root, but 
gives complex roots only for Figure 9-Al, where ? = 0.7 and > 50. 
Each of these figures shows a very fast initial transient, there 
is no apparent contribution from the small real root (Table I, 
column y ) t which may have been cancelled by the zero of the 
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transfer function. Comparison of these four initial transients 
shows very good agreement with the root values of Table I for 
low leakage (Figure 9-Al) . The initial transient shows a small 
overshoot, about as expected for roots with c = 0.7. 
increases (Figures 9-A2, A3, A4) the transient does not overshoot 
(roots are all real) and the settling time increases which 
agrees with the change in root values in column S of Table I. 

Continuing the analysis of Figure 9, the z-displacement and 
plenum pressure curves show an approximately exponential 
variation which is very slow but speeds up with increased leakage. 
This agrees qualitatively with the increase in value of the 
smallest real root (column y, Table I). The pitch angle variation 
of Figure 9 shows a substantial variation, but the short 
duration does not supply enough information for further analysis. 
We therefore proceed to Figure 10. 

Let us start the discussion of Figure 10 with the pitch 
angle response curves, 10-Dl, D2 , D3. On all of these curves it 
is clear that the pitch transient is an exponentially damped 
sinusoid plus another exponential response. The frequency of 
this pitch mode is approximately 5 rad/sec, and increases slightly 
with increased . We cannot correlate this with the linear 
model, for it does not include the pitch mode. The second 
exponential component, however, appears to be due to the small 
real root of the heave mode (Table I, column y) and its time 
constant seems to decrease with increasing K^, which agrees with 
the root variation in Table I. 
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From Figure 10 it is also clear that the pitch mode couples 
into c.g. acceleration, z-displacement , and plenum pressure. 

This coupling increases with K^. The small real root is also 
seen to contribute — it is obviously dominant in the z-displacement 
and pressure curves, and can be recognized in the c.g. acceleration 
curve of Figure 10-A3. 
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FIGURE 9-Dl. Step Weight Transient 

= 6.28 1 sec, 
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FIGURE 9-C4 • Step Weight Transient 

= 21.99 t^ < 1.0 sec. 
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= 14.66, t^ = 5.0 sec. 
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FIGURE 10-C3. Step Weight Transient 

= 21.99, t^ = 5.0 sec. 



0C3 










■j- ^ 



UER3U3 TIME 



o 




V. CONCLUSIONS 



A few specific conclusions may be drawn from the 
results of this study: 

1, The major effect of increasing air flow rates is 

to attenuate high frequency components in the response 
of the craft. This is clearly shown by the high 
frequency response of both the linear and 6 D.O.F. 
models. It is also shown by the step weight removal 
transients. Therefore, to reduce the vibrations due 
to vertical accelerations at high encounter frequencies, 
the air flow rate should be maintained as high as 
practical , 

2 , The pitch mode couples strongly into the vertical 
motions over the encounter frequency range of 1 to 10 
rad/sec, as shown by the 6 D,0,F, model. Therefore, 
pitch motion should be added to the linear model to 
account for the changes in air flow leakage rate 

due to the variations in the pitch angle, 

3, The linear model (Heave motion only) is sufficiently 
accurate for many studies and should be used where 
applicable to aid in the analysis and design of CAB 
craft characteristic, 

4, The linear model reduces both the CPU execution time 
and the operating class in the job queue resulting in a 
reduce interaction (turnaround) time. 



66 



5. The complex sea method of developing the frequency- 
response of the 6 DOF model using the Fast Fourier 
Transform provides approximate results that can be 
quite inaccurate unless special care is exercised 
in its application. 
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VI. RECOMMENDATIONS 

One of the long range goals of the analytical studies 
of the CAB craft at NPS is to develop a simplified model 
that will provide a reasonable computer execution time for 
control system design and evaluation in the complex sea 
operating environment. We feel that we have a reasonable 
start in this direction and have the following recommendation 
for future studies: 

1. The pitch mode should be modelled and incorporated 
in the Linear Heave-Only Model. Mr. McIntyre in 
Reference 7 has made the first effort in this 
direction and a current study by Mr. Barnes (Reference 8) 
has improved on the model and shows promise of 

r providing the additional information missing from 

the Heave-Only Model . 

2. Upon the completion of the Linear Pitch-Heave Model 
and its regular sea frequency response validation 
with the 6 D.O.F. model, the complex sea frequency 
response should be computed and compared for both 
accuracy and speed of computation. 

3. Continue experimentation with both software and hardware 
techniques for the analysis of the complex sea results 
both for the 6 D.O.F. model and the actual XR-3 craft 

to provide improved results of the complex sea frequency 
response analysis. 
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4. Develop deterministic input tests for the XR-3 
craft for validation of the SESCOMP model. The 
weight removal transient is one such test that 
shows promise for validation of the Pitch-Heave 
mode equations. 

5. Investigate effect of the addition of a membrane 
in the plenxom chamber to reduce the variations in 
the vertical motions of the XR-3 craft. The 
actual XR-3 craft as well as the 6 D.O. F. program 
and the linear pitch-heave model would be used 

in this study. 
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APPENDIX A. DESCRIPTION OF THE SIMPLIFIED XR-3 MODEL 



A. Nonlinear System 

In this model of the XR-3, pitch variations are reduced 
to zero by assuming 1) the center of pressure (C.P.) directly 
under the center of gravity (C.G.) , 2) the sidewall symmetrical 
about the C.G. with uniform rectangular cross-section from 
bow to stern and 3) the pitch moments of the seals and 
aerodynamics cancel each other. The lift force of plenum 
pressure and sidewall buoyancy oppose the craft weight. All 
other lift and drag forces are neglected since constant speed 
conditions are assumed. 

In addition it will be assumed that the rear seal maintains 
a constant leakage area, that is, the seal follows the water 
level at a fixed separation. 

Figure A1 shows a top and sideview of the simulated CAB 
craft. All dimensions are chosen to approximate those of 
the XR-3 craft in order to compare the results of this 
analysis to the L & M program solution under similar operating 
conditions. 

The equations of motion for this system are given below. 



Orifice Leakage Rate, *3^^^ 



C A„ 
n l 




cu. ft. 
sec. 



( 1 ) 



Fan Map Input Rate, 



n 



[Oi (o) 



k Pw 1 

q b J 



cu. ft . 
sec . 



( 2 ) 
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FIGURE Al.-Top and Side View of Simplified XR 



Absolute Plenum Pressure, P, = P^ 
(Adiabatic Process) “ 






Psf 



(3) 



Plenum Volume, Vj^ = (V^ - cu.ft. 



lbs , 



Plenum Air Flow Rate, Mj^= ^out) ii^ 



(4) 

(5) 



\ M M 


M / 


sec . ^ 


■'6) 


Heave Velocity, Z = 


ft. 
sec . 




(7) 


Heave, Z = /z dt + 


Z (0)1 


ft. 


(8) 


L*^o 


J 






Plenum Pressure Lift Force, F 

P 


= "^b^b 


lbs . 


(9) 


Buoyancy Lift Force, 2 (A 




] lbs. 


(10) 


Plenum Gage Pressure, Pj^ = (Pj^ “ 




Psf. 


(11) 


Draft, = Z + 


Zs 


feet 


(12) 



The system parameters and constants are listed below. 



Adiabatic process coefficient. 
Leakage area. 

Leakage orifice coefficient. 
Air density. 

Atmospheric pressure, 



Y = 1.4 

C = 0.90 
n 

p = .002378 
a 

P = 2116. 
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Plenum area. 

Empty plenum volume, 

Craft weight. 

Keel line area of sidewall. 
Draft, initial condition. 
Density of water. 

Gravity acceleration. 

Center gravity location. 
Number of fans. 

Steady state fan output. 

Fan map slope. 



= 200 . 

Vn ^ 383. 

W = 6007 
A = 75/4 
£^(o) = 0.5 

p = 2 . 

g = 32. 

Zs = 2.5 
n = 5 . 



Q^(o) = 35. 
kq = 1.0 



The schematic diagram of signal flow is shown in Figure A2 . 
The coordinate system and initial conditions for heave and 
draft are shown in Figure A3. 




Figure A3 - Coordinate System and Initial Condition 
of Heave and Draft. 
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Z (0) 



Z(0) 




FIGURE A2. — Schematic Diagram of Signal Flow. 
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It should be noted that the range of heave motion is 
-2.5 < Z < -0.5 since -2.5 feet would put the water level at 
the keel line and -0.5 feet would put the water level at the 
plenum roof. Draft is measured from the keel and would have 
the range 0 < < 2.0. Also note that if the craft "drops" 

into the water, increases and Z becomes less negative, that 
is, Z increases in the positive (downward) direction. 

At steady-state condition, the plenum pressure assumed 
is - 24. PSF producing a ^^( 0 ) = 4800 lbs. For the 

draft of = 0*5 feet, the buoyancy force = 1200 lbs. 

The air mass in the plenum for the above steady-state 

where q. (o) = q . (o) = 55 cu. ft. /sec. is M, (o) = 0.71878 
^in ^out b 

pounds . 

B. Linear System 

The linear system equations are arrived at by 
application of the Taylor Series expansion about the 
steady-state operating point. The development is shown below 
for the two differential equations. 



a) Steady State and incremental variations are noted as follows 
using heave as an example. 

2 = 2 (o) + 62 
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Equation (1) , 



let W/M - W„; F^/M 
m p 

and F^/M = F^^ 



F 



pm 



then, 

2(0) + 62 = W/M - (Fp^(o) + 6Fp^) - (F^^(o) + 6F^^) 



or. 



• • 






(la) 



Equation (2) , let 



then , 

V^(o) + 6V^ = nQ. - nk (Pj-(O) + 6P ) 

m m 4 o o 



C A 

n £ {— ) 



2P. (o) . 

C A„ is (—2 ) ^ 

n c p 



{— } 5P, 



or. 



• ^ a. h — 

6V = -(nk^)5P, ^ 5P 

^ ^ “"a 2P^(o) ^ 

D 



(2a) 



Next, the auxiliary algebraic equations are linearized. 

F^(o) + 6F^ = A. (P. (o) + 6P, ) 

p p b b b 




(3a) 
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(o) + 6F^ = 2A^ p^g U^(o) + 6^^) 



»w9> 



(4a) 



Eg. (5a) and (6a) have already been developed they are: 



and 



6q. = -(nk ) 6P>^ 

^in q b 


(5a) ; 

1 

1 




C A p p 1 

'4out 0^ b 


1* 

(6a, 1 
j 
1 



For equation 7a 



viV . 



(V (o) ) 1 

1 1 6V_ _ m 



V^(o)-] 



■ (V^(o) ]2 



6P, = P Y ( ) 

b a 



,-i 



m m — m b 



V, (o) V^(o) V, (o) V, (o) 

b m b b 




and equation 8, 9, and 10 follow. 
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\(o) + 6V 



b 



V 

n 



Ab (£^(o) + Sl^) 





= - 


(8a) 


(o) + 


= ^b<°> + «'’b - '■a 




'^b 


= 


(9a) 


(o) + 


St. = z(o) + 6z + Z, 

U O 




d 


- 5z 


(10a) 



Combining (la), (3a), (4a), (7a) and (8a) as follows: 



• • 

6z = 





6 z 



recall that = 5z 

a 



6z 





Pb(o) 



1 6V^ - 1 

V^(o) ^ V (o) 

m b 
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6 z = - 



6z + 

/ \ MV. (o) /- 

D 



YPb(o)Ab| 

. MVj^(o) J 



m 



Also for equation (2a) 



6V = 
m 



nk + / ^a \ ^ yPb^°^ 

»^a V2Fb(o)/ -* 



6V 



m 



6V. 



_V^(o) V, (o) 
m D 



6V = 
m 



^ TPb<o)[( 

Pa ^ 2 P (o)^ ^ 

D 



+ / ^ \ 62 



and for equation (7a) 



6Pb = YPb(°) 



• 1 6V + b 
m 



6 z 



m 



Vb(o) 



Now define the state vector 



X A [6z 6z 6V ] 
— = m 



and output variable, y = 6P. with 



y A Cl 1 Xi + Cl 3 X 3 



Then the state equations are 
xi = X2 

X 2 =321 5Ci + a23 X 3 
^3 = ^31 ^1 + ^33 ^3 
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, yPi^(o)a, j yPi_ (o) 

where Cii = b b and C 13 = ' b 

v^(o) v^(o) 

/2A ,p g , yPw(o)A. 
a?i = -[ 4 w^ + ' b b J 

V M MVj^(o) / 






33 1 = - yP.- ( 0 ) nk_ + ^n^£ 


"a 


J 

1 


D Cl 


J 2 Pj^(o) _ 


a 33 = - YPh(o) nk + ^n^-^ 


"a 


L "a \| 


2?b<0)J 


V„(o) 

m 



The characteristic equation for this system is given by 

#► 

(4 1 - A) where 
A = 



o 
a 2 1 



1 

o 



o 

^2 3 



^3 1 



Then (41 - A) = 



o a 3 3 

4 -1 o 

-321 ^ ~^23 

—a 3 1 o 4 — a 



33 



and the characteristic equation^ is 
3 2 

4 - 333 4 - 321 ^ + (^21 ^33 " ^23 ^3 1 ) = O 



Using the given numerical values the characteristic equation is 
3 2 

4 + 62,4 + 21834 + 789 = o 

This cubic equation will factor into one real root near the 
origin and a complex pair far removed from the jw axis. (See Sec. 
IVA) Thus a step weight response will have a complex pair of poles 
which can be approximated by the quadratic form. 
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2 

-6 - a33'i~^21“^ 

For the given numerical values it can be shown that a 
rough approximation for the coefficients are: 



-a 33 = nk^^ T P^(o)/V^(o) 



eq 



-a 2 i = Y P. (o)/V (o) 



M 



Using the above approximation yields a natural frequency 



of : 



Wn - - 



^ Y Pj^(o)/Vj^(o) 
M 



w = K , 
n n ^ 



Pb(o)/Vb(o) 



(13) 



where \j y/M 



and a damping factor 



^ = ~^33 



2w 



n 



n P^^ (o) / V^J^^(o) 

4 ^b^ Y 

M 



C = K, 



\| 



Pj^(o) Vj^(o) 

V2 (o) 
ra 



(14) 



where 



= 



nk 

c 

2A, 



\l Yin 



82 



APPENDIX B. COMPUTATION OF THE WAVE INPUT TRANSFER FUNCTION 



The block diagram of Figure 1 given in Section IIIA can be 
drawn as a signal flow graph as shown in Figure B-1. To obtain 
the system transfer functions, the Mason Gain Rule (MGR) can be 
applied as follows: 

A = 1 — (Lj + L2 + L3) + L2L3 + L2L4 

where 



Li = -k2 


k 3 k 4 


CM 

1 

00 

Ai 


ro 

1 

II 

<M 


kg k7 




L3 = -ks 


CM 

1 




L4 = -kg 







For the incremental heave (draft), 6z, the MGR yields 
for the wave input, 6W 

6z = I (Pi Ai + P 2 A 2 ) 6W 

where 

Pi = “kflA ^ Ai = 1 — L 2 

P 2 ~ “k 1 k2 ks k If 4 ^ A 2 “ 1 

thus , 

6z = i [(-ko4"2) (1 + kjksky^"*) + (-kik2k3k44"^) (1)] 6W 

where the gain coefficients for any particular operating 
conditions are: 



ko = 



2W V 
s w 

cj • 



sin 
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1 



I 

( 



m 







FIG B-l SIGNAL FLOW GRAPH OF HEAVE EQUATIONS FOR WAVE INPUT 



^2W V , 
p wi 







\“i / 


ka 


= 


(O) 


kg 


= 


YPb (o) 


k4 




Afa/M 


ks 




2A^p^g/M 


ks 


= 


V <3 1 



sin 



2_i 

pa 





k7 = 1/V^ (o) 



k3 



kg 



I 8M 



By applying substitution and reduction, the transfer function 
becomes , 



6z _ -(a-6 + b) 

<Sw + C4^ + d4 + e 



with 

a = (ko + ki kz kg k 4 > 

b = ko kg ks k 7 

c = (kg ks k? + kg) 

d = (ks + ka kg k4 ka + kg ks k? kg) 

e = kg ks ks k7 
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It should be noted that the transfer function of draft 



to wave height has a numerator whose coefficients "a" and "b" 
are functions of the incident wave, co^, and therefore will 
change in value for each frequency of the sinusoidal sea input. 

In addition, it should be noted that the incremental heave 
(draft) is 180 degrees out of phase with the wave input as 
would be expected since a positive value of wave height produces 
(1) increased buoyancy and (2) reduction in plenum volume which 
increases the plenum pressure and thus the lift force. 

Further, it is observed that the denominator polynomial 
is of third order as would be expected since the characteristic 
roots of this system are obtained by solution of the equation 

4^ + c4^+d4 + e = 0 

as given in Section IV-A. The transfer function of c.g. 
acceleration, 6z, to wave height, 5w, is obtained from the 
signal flow graph as follows. 



5z 

6W 



-kik2k3k4(l) - ko(l + k3ksk74~M 



A 



6z 

6W 



-4 ^ (a4 + b) 
4^+c4^+d4 + e 



where the coefficient a, b, c, d and e are the same as those in 
the draft-wave height transfer function. 
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The minus sign appearing in this transfer function represents 
the 180 phase shift which is also shown in the draft-wave height 
transfer function as would be expected. 

In like manner, the transfer function of plenum pressure, 

6Pj^, to wave height, 6W is 

kikakad + ks^"^ + kg^”^) - kokakaks^”^ 



_ 4{a^4^ + b^4 + cM _ + bM 

A A 

= ki ka ka 

= ki ka ka kg 

- ka ka (kiks - kokg) - zero 



6W 



where 



a* 

b^ 
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APPENDIX C. RELATIONSHIP OF ENCOUNTER 



FREQUENCY TO INCIDENT WAVE FREQUENCY 
FOR THE XR-3 CRAFT 



Booth (in Appendix A of Reference 4) gives a brief 
development of the equation for encounter frequency and 
for the case of ahead seas, the incident frequency, is 
shown to be 











- 


•H 

3 


= 


g 

2V 




— 1 

* — 1 

(U 

3 

+ 

» — 1 



This relationship was plotted by Booth and is repeated 
below in Figure C-1. For the XR-3 craft operating at a 
steady-state speed of 30 knots, V = 50.67 ft/sec., the 
incident wave frequency becomes 

0)^ = 0.318 [jl + (6.3)ojg - l] 

An approximation to this functional relationship that is 
useful over the range 

1 < 0 )^ < 20 rps 

is 0 )^^ = Wg/2 for the 30 knot speed fo the XR-3 craft. 
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Figure C-1. Wave Frequency vs. Encounter 
Frequency 
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